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Abstract 



Collective motion is abundant in nature, producing a vast amount of phenomena 
which have been studied in recent years, including the landing of flocks of birds. 
We investigate the collective decision making scenario where a flock of birds de- 
cides the optimal time of landing in the absence of a global leader. We introduce a 
simple phenomenological model in the spirit of the statistical mechanics-based self- 
propelled particles (SPP-s) approach to interpret this process. We expect that our 
model is applicable to a larger class of spatiotemporal decision making situations 
than just the landing of flocks (which process is used as a paradigmatic case). In the 
model birds are only influenced by observable variables, like position and velocity. 
Heterogeneity is introduced in the flock in terms of a depletion time after which a 
bird feels increasing bias to move towards the ground. Our model demonstrates a 
possible mechanism by which animals in a large group can arrive at an egalitarian 
decision about the time of switching from one activity to another in the absence of 
a leader. In particular, we show the existence of a paradoxical effect where noise 
enhances the coherence of the landing process. 
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1 Introduction 



Scenarios involving collective decision making appear to be ubiquitous in sev- 
eral fields including animal behaviour [1] and social sciences [2j . Several stud- 
ies [311^ have investigated the cases where any group unanimously decides to 
choose one of the many available options. It is generally expected that there 
will be diflFerences in the motivations of the members at the time when mak- 
ing choices. However, in spite of the differences consensus is seen to evolve. 
Different mechanisms have been suggested to account for this type of process 
[51I61I7] . The related phenomenon of collective motion in animal groups has also 
been extensively modelled [HllQlfTO] by considering animals as point-like parti- 
cles (so-called self-propelled particles). In these models it is assumed that each 
individual moves with constant speed while tending to align with its imme- 
diate neighbours, for low level noise, giving rise to a globally ordered state 
replicating the motion of fiocks where animals move in the same direction. 

Our object of study is the phenomenon during which the animals moving in 
groups seemingly make unanimous decisions on the time and choice of per- 
forming activities plfTT] even in the absence of global leaders. Examples include 
takeoff of swarm of honeybees from nest sites [12j, activity synchronisation in 
Merino sheep [13j, collective movement of white-faced capuchins [Hj, group 
departures of domestic geese [T5] and departure of Argentine ants from feed- 
ing site |16j. Such a paradigm could be generalised to human behaviour as 
well, e.g. where to stop to rest, when making an excursion with a group. The 
important extra feature of taking into account group motion during collective 
decision making is that in such an approach the neighbourhood with which a 
consensus is to be achieved is dynamically changing in a realistic manner. 

In this report we model the process of landing of bird fiocks performing forag- 
ing fiights as a typical example of collective decision making. We regard the 
birds as self-propelled particles, the only difference between the birds being an 
a priori value corresponding to the heterogeneity in motivations. Throughout 
this paper we shall use the terms "birds" and "particles" interchangeably. In 
a recent model for collective landing [TT] the birds are assumed to move under 
the action of different social forces [18j. In addition, the internal state of each 
bird is characterised by a continuous variable called landing intent such that 
the internal state of each bird is directly coupled to the internal state of its 
neighbours. Another model [19j allows the motivation of individual birds to be 
infiuenced by only the observable variables of their neighbours, like velocity 
and position. This assumption simplifies the description of the landing pro- 
cess, since it does not presume very sophisticated channels of communication 
during fiight. The model discussed in this report has also been formulated 
along similar lines. The decision of a bird is infiuenced by the state of mo- 
tion of its neighbours and an a priori intent (a pre-assigned constant) taking 



2 



into account its endurance that acts through a noise variable. This makes the 
landing a stochastic process even if we consider a single bird only. 

Our aim with this report is to quantitatively show that a group without leaders 
can perform a synchronised landing following very simple rules, even when a 
single bird does not have information about the whole of the flock. One of the 
important new features of our approach is the possibility of tuning the level 
of perturbations that the particles are subject to during the landing process. 
This approach also allowed us to uncover an interesting effect. According to 
our simple model, there is an optimal level of perturbations resulting in the 
smallest spread of the landing times of the individual birds, meaning that in 
certain cases the noise can enhance the coherence of the landing process. 



2 Model 

2. 1 Basics 

In this section, we will describe the basic concepts of our model, leaving the 
majority of the technical details to the next section. We formulate the dy- 
namics of a bird flock along the lines of social forces. These forces are similar, 
but not the equivalent of forces considered in earlier studies Pl20ll21j . In the 
model, we consider a flock that performs a horizontal flight about an average 
height until it decides to land. This allows us to handle the horizontal and 
vertical motion of the birds separately. 

Let us flrst regard the horizontal motion. Setting z = as the only natural 
boundary, representing the ground, the horizontal flight is performed parallel 
to the xy-plane. There are four social forces acting upon the ith bird horizon- 
tally, (i) f^y - an averaging on the velocities of other birds in its neighbourhood, 
(ii) f^y - an evenly distributed noise i.e., a random perturbation on its veloc- 
ity, (iii) fjy - a repulsive force, and (iv) f^y - a cohesive force. The force f^y is 
responsible for aligning the velocities of neighbouring birds in the flock. For 
this purpose we deflne the neighbourhood {Ni) of the ith bird as those birds 
that are encompassed by a cylinder of inflnite height and radius R centred 
on the ith bird. While in flight a bird would like to maintain some separation 
with other birds. This is accounted by the repulsive force. This force is char- 
acterised by a certain radius of interaction rf, which can also be considered 
as the eflFective "size of a bird" and is much lower than the interaction radius 
characterising Mi. The purpose of the cohesive force is to keep the flock from 
breaking apart. It is described by another inflnite cylinder with a radius D 
which is much larger than that of A/^, and is centred on the centre of mass 
(CoM) of all birds. If a bird strays out of this cylinder it is propelled towards 
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the centre of mass. 



The vertical motion also consists of four social forces: (i) - an averaging on 
the vertical velocities of the neighbourhood Afi of the ith bird, (ii) - a force 
that tries to keep the birds about a given height, (iii) - a repulsive force 
like in the vertical direction, and (iv) - a noise. The force is such, that 
around a given height h there is a regime of width Ah where in all practical 
sense, a bird moves freely, but outside of that, the bird is repelled toward h. 
The most important part of the model is the vertical noise, since this is the 
force that facilitates the landing of the birds. This is done as follows: each 
bird is assigned an a priori value, a U depletion time, which represents the 
time at which the bird starts to feel the depletion of its energy reserves, and 
would increasingly wish to land. In general these depletion times will depend 
on several external or internal conditions such as the energy reserves of the 
birds [22], stamina, health, willingness to fly and other things and thus will be 
in general different for different birds. As such, these values are drawn from a 
Gaussian distribution of a mean /x and a variance a^. When a bird i reaches its 
ti^ the originally evenly distributed noise starts to get increasingly biased 
to facilitate landing with a characteristic time r. Landing occurs when this 
force and the averaging force overwhelms f^^ allowing the birds to land. Thus 
the equations of motion of the ith bird in the horizontal direction are: 



pum ^ fa ^ p ^ fc ^ p 

xy,^ xy,2 ' xy,^ ' xy,z ' xy,z \ / 



and 



pum 

r,y,,(t + At) = r^y^t) + Vj^At, (2) 

\ xy,i I 



where r^y^i is the position of the ith bird in the xy-plane, and is a constant, 
so that the birds move with a velocity of constant magnitude. Very similar 
equations hold for the vertical motion: 

pum ^ p i I fh I p . (q.) 

J z.i J z.i ' J z,i ' J z,i ' J z.i \ J 



^sum 

z,it + At) = z,it) + v^At, (4) 

\Jz,i I 

where Zi is the vertical position of the bird. As seen in the next subsection, the 
magnitude of these forces are parameters of the model, which seemingly are 
a lot in number, but only some of them are actually relevant to the problem 
of landing. Therefore we study the behaviour of the system as a function of a 
few parameters, viz. the variance of depletion times a, the number of birds A^, 



and the coefficients of f^- and f^- (described in 2.2). The characteristic time 
T is only important in the sense that it must be of the order of magnitude of a 
for non-trivial behaviour of the model. In addition, all the differences between 
the birds are contained in a. 
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2.2 Details 



In this section we provide the exact mathematical expressions for the forces, 
and other details concerning the model. Noting that Vxy and Vz stand for 
the instantaneous velocities in the horizontal plane and the vertical direction 
respectively, the averaging forces are given by: 

fxy. = (Vxy)x (5) 
= (6) 

where ( is the average over the birds in the neighbourhood Afi of the ith 
bird. We also assume that the averaging includes birds which have already 
landed (in jVi) so that their influence is taken into account. The eflFect of the 
landed birds on the averaging forces would be the same as that of those moving 
vertically downwards. The repulsive forces follow the same scheme in all three 
dimensions. The explicit form for the z-direction is as follows 

N 

f z,i ^ J f z.ii l'^) 



where 




, -- if < \zi — Zn\ < d 



otherwise. 



Here A is the strength of the repulsive force. The cohesive force on the ith 
bird is given by the following equation: 



n -f l^ACoMl ^ D 

^ iJ- |lxy,i I ^ 2 

-B {\r^yf^\ - f ) otherwise. 



(9) 



where r^t^f^ is its distance from the CoM and B is the strength of the force. 
The form of the force ensures that if a bird is farther from the C'oM than D/2 
it feels an attraction towards the CoM. The social force describing the will to 
stay at a given height h is given by 



20 



1 + 



tanh|^(|z.-/.|-^)| 



sign(z^ - /i), (10) 



where C is the strength of the force. The attraction f^- towards the preferred 
altitude h is actually weak within a region of width Ah and is strong outside 



(see flgure la). 



We now deflne the nature of the noise in the vertical and the horizontal di- 
rections. At any instant of time the ith bird is influenced by a horizontal 
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Fig. 1. (a) The force against z. The plot shows the small forceless regime around 
h and the fast strengthening of the force outside of that, quickly saturating to a 
constant, (b) The plot shows the time evolution of the probability distribution of 
f^. Equation [13] is the explicit formula for generating values of this force. 



noise: 



fxy.W=/3exy.(t), (11) 

where (^xy,z is a unit vector on the x^-plane whose orientation is taken to be 
random and /3 is the strength of the noise. For times less than ti the vertical 
noise /z i is given by: 

f:At) = ^^.At), (12) 

where Cz,i{t) is randomly chosen from a uniform distribution on the interval 
[—1,1] and a is the amplitude of the vertical noise. After time t crosses ti the 
form of gets modified to: 



a — a J I + 4tan6' + ^tan^O - 8^z,i(^)tan6' 



2tan6' 



(13) 



where 



tan^ 



1 + exp 



t-U 



The above form ensures that after ti the noise becomes continually biased 



towards the downward direction (see figure lb). The characteristic time scale 
for this biasing is given by r. 

For setting the units of measurement we choose i? = 1 and At = 1. The 
summary of the parameters, and their typical values, as used in our simulations 
are provided in table [T] As initial condition the birds are assigned randomly 
oriented velocities and are uniformly distributed in a cylinder of diameter D 
and of height A/i. 

For a rough fitting of the model to real fiocks we took data from [23j . Based 
upon their measurements it seems reasonable to take 10 m/s and i? 2 m. 
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Table 1 

The parameters with their brief explanation and values as we most often used them 
in our simulations. We find that the most relevant parameters to the problem of 
landing are a, (7, a and N. The values of the parameters were chosen so to obtain 
a biologically relevant landing scenario. For setting the units of measurement we 
choose R = 1 and At = 1, which makes the horizontal interaction radius (the 
radius of Af) 1. 



Par am. 


Description 


Value 


Dimension 


a 


the coefficient of the vertical noise 


0.2095 


R/At 


C 


the coefficient of 


0.81 


R/At 


a 


the standard deviation of the distribution of ti 


5000 


At 


N 


the number of birds 


300 


- 


V 


the speed of the birds 


0.1 


R/At 


D 


the horizontal diameter of the flock 


20 


R 


d 


the "size" of the bird 


1/6 


R 


h 


the optimal height of flight 


100 


R 


Ah 


the effective width of the flock 


2 


R 


T 


the timescale of energy depletion 


5000 


I /At 




the mean of the distribution of U 


50000 


At 


A 


the coefficient of and f^y 


100 


1/At 


B 


the coefficient of f^y 


100/3 


I /At 


P 


the coefficient of the horizontal noise 


0.05 


R/At 



which makes At ^ 0.02 s, since ^' = 0.1 i?/At in our model. As a rough 
estimate from the simulations, we can take the time needed for the flock to land 
as 5000 At. This means that our simulated flocks flying at the altitude of 200 m 
land in about a 100 s. Taking into account, that our model is very simplified 
compared to the actual dynamics of a bird flock, the different parameters of 
the model are in a reasonable range for real-life bird flocks. 



2.3 Quantitative characterisation of the landing 



To quantitatively characterise the cohesion of the landing, we introduce a few 
quantities. We measure the spatial cohesion of the flock by taking the standard 
deviation (cTxy) of the horizontal coordinates from the centre of mass as the 
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^ 2 

E (rxy. - r^r) (14) 



N ■ 1 



following: 

^xy = 

where 

CoM ^j=i ^-^yd 
xy ^ 

We sample this quantity at two different times, once during flight (cr^y) and 
once after the whole flock has landed (cTxy), so that their ratio forms a measure 
to characterises the degree of coherence of the landing in space. 

To characterise the temporal coherence of the landing we measured the fol- 
lowing two quantities viz. the standard deviation of the times at which the 
different birds have landed (ctl) and the time (Tqq) that elapses between the 
landing of the 20% of the flock and the landing of the 80% of the flock. We 
additionally calculated the latter quantity in the case when the coupling be- 
tween the birds was set to zero, i.e., when = denoted by Tqq. To obtain 
normalised measures characterising the temporal coherence we took ctl/ct and 



2.4 Ordering horizontally 



For obtaining a biologically relevant scenario, we adjusted the coefficient of 
the horizontal noise /5 to an appropriate value. As can be seen in figure [2| 
that starting from the totally coherent, no noise situation, the increase in /5 
gradually destroys this coherence leading to a situation, where the centre of 
mass of the fiock does not move at all. In this case we see the boundaries of 
f^y if we plot the tracks of the birds projected unto the xy-plane. We chose 
to investigate the regime where the horizontal tracks would be similar to the 



one illustrated in figure 2b 



2.5 Assay of the landing 



After choosing the parameters for horizontal motion we observe the time evo- 
lution of the z-coordinates of a similar fiock in figure |3| The two states of 
fiight, horizontal motion of the fiock at a height around h and the landing are 
clearly visible. To see whether the landings are synchronised or not we have 
plotted the percentage of the landed birds in the fiock as a function of time, 
in the presence of coupling, in the absence of coupling {f^ = 0) and in a mean 
field case (meaning the radius of Af is infinity, not R) in figure [ij 

As we can see, the transition in the coupled case is much sharper than in the 
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20 40 60 80 



(a) /? = 




-50 -40 -30 -20 -10 



(b) /3 = 0.03 



-20 -15 -10 -5 5 10 15 20 



(c) /3 = 0.05 



-20 -15 -10 -5 5 10 15 20 



(d) /3 = 0.07 



Fig. 2. The tracks of the birds projected onto the x^-plane with different levels of 
horizontal noise. The different colours represent different birds, where 10% of the 
birds were plotted with longer tracks and the rest of the birds with just the end of 
their paths. From (a) to (d) we observe a transitions from an unrealistic straight line 
of movement, through the biologically relevant phase to the total loss of coherent 
movement. 
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Fig. 3. The z-coordinates of the different birds against time. The flight around 
the height h gives way to collective landing over time as more and more birds are 
increasingly biased towards moving downward. The parameters were chosen as of 
table [H 

uncoupled case, showing considerable synchronisation among the birds. It is 
also notable, that increasing the radius of interaction to infinity does not make 
the landing process relevantly sharper, it merely decreases the time needed to 
make the decision to land. Coherent landing arises from the interplay of three 
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time 



Fig. 4. The percentage of landed birds as a function of time. The red curve corre- 
sponds to the case when couphng between the birds is absent, i.e., = 0, the green 
one corresponds to the coupled case, while the blue curve is the mean field case, i.e., 
where the radius of Af is infinity instead of R. It is clearly seen that in the presence 
of coupling, the landing is much sharper viz. the synchronisation among the birds is 
much greater. It is also notable, that increasing the radius of interaction to infinity 
does not make the landing process relevantly sharper, it merely decreases the time 
needed to make the decision to land. The parameters were chosen as of table [l] 

Table 2 

The normalised quantities describing the coherence in space and time, with the 
parameter values from table [l] Tqq is the value of Tqq for the zero coupling case and 
is the values of a^y during flight. 



Quantity 


Value 




0.387 ±0.007 




0.288 ± 0.007 




0.639 ±0.017 



forces: the vertical noise and the averaging force. The magnitude of the 
bias in the noise increases with time to ultimately overcome f^^ but due to 
the averaging force the individual biases are, in a sense, averaged over the 
neighbouring birds. Thus for an individual bird it becomes "harder" to land 
when the bulk of the flock would still want to fly, and becomes "easier" when 
the latter wants to land. 



In table 2 we have summarised the quantitative analysis of the landing. For 
this we measured the quantities introduced in section |2.3| in around hundred 
independent runs of the simulation, and averaged over the obtained data. Note 
that the normalised quantities describing temporal coherence show consider- 
able synchronisation in the flock. We stress that this happened while a single 
bird had information about only a few neighbours, and not the whole flock. 
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Fig. 5. The values of the most important quantities describing the collective landing 
for different number of birds, while keeping the density of the birds constant with 
the changing of diameter D (all other parameters are as of table 1). Note that 
Tqq/Tqq is well below unity and this property persists for much larger flocks than 
the one used in most of our measurements. 

In figures [5| |6] and [7] we plot the quantities mentioned in table [2} as functions 
of the number of birds in the flock, the standard deviation of the distribution 
of and the magnitude of the vertical noise, respectively. While changing 
the number of birds in the flock, we keep the density of the birds constant by 
appropriately changing the horizontal diameter D of the flock. Note that our 
previous statements about temporal coherence hold for considerable changes 
of these attributes of the flock. 

In figure [5] we see that decreasing the number of birds, while keeping the ratio 
of the number of birds and the volume set by the boundaries of the forces 
and a constant, makes the flock less coherent. This is due to the fact that 
with the decrease in the number of birds, the fluctuations in the number of 
neighbours, with respect to the individual birds, become important. In figure|6] 
we see that for small as the value of ctl/ct grows well above unity. This 
is because there is an inherent difference in landing times, due to the birds 
starting from different heights and thus arriving at the ground at different 
times, and also due to the fact that repulsion prevents the birds from landing 
on top of each other. The effect of this inherent difference becomes noticeable 
compared to low a values. 

Figure [t] shows that for any given coefficient of and (7, there is maximum 
of temporal coherence or, alternatively a minima of Tqq/Tqq and ctl/ct, as a 
function of the noise. The interplay between and the vertical noise be- 
comes apparent from figure [7b| The location of the minima as a function of 
a (the coefficient of the vertical noise) shifts to right, and its value increases 
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Fig. 6. The values of the most important quantities describing the collective landing 
for different standard deviation of the depletion times (all other parameters are as 
of table 1). Note that Tqq/Tqq is well below unity. For small values of a, the ratio 
(Jl/ct grows considerably due to the fact, that factors such as the size of the birds 
and the height of the flock, become important in deciding the time taken by the 
whole flock to land. 



(in (Jl/ct) with the increase of the coefficient of f^. To explain this, let us 
increase a from zero as we keep C a constant. At small values of a we see 
that there is no landing, as the noise is not strong enough to overcome f^. 
When a is larger (right side of figure 7b) the temporal coherence of the fiock 
decreases with the increase of a, as one would intuitively think. Between these 
two there is regime with a non-trivial minimum of ctl/ct (maximum of tem- 
poral coherence), where the increase in the magnitude of the noise actually 
increases temporal coherence within the fiock (see e.g. noise-induced ordering 
[24j . stochastic resonance [25]). The increase of C decreases the maximum of 
temporal coherence the fiock can reach and increases the magnitude of noise 
needed to reach this maximum. 



Our model, naturally, uses assumptions which are results of simplifications 
of the true complexity of a fiock of birds. We assume local interactions (that 
these dominate the landing process) in addition to the a priori assigned intent 
of the birds. A possibility would be to include the reaction of a single bird 
to the behaviour of the whole fiock (taking this into account would need a 
number of further arbitrary assumptions), but we have restricted the model 
to mostly local interactions, with one bird interacting with an estimated 8-12 
birds on average. 
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Fig. 7. (a) The values of the most important quantities describing the collective 
landing for different magnitudes of the vertical noise (all other parameters are as of 
table[T]). For values of a slightly less than the smallest shown on the graph no landing 
occurs, (b) The temporal coherence of the flock as a function of the magnitude of 
the vertical noise and (all other parameters are as of table [l]). The place of the 
minimum in a shifts to the right, and the value of the minimum increases as 
becomes stronger. For values of a slightly less than the smallest shown on the graph 
no landing occurs. 



3 Conclusion 



In this report, we have investigated a model where the concepts of collective 
decision making and collective motion are intertwined. In particular, we have 
introduced a simple phenomenological model for the landing of a flock of birds. 
In the model birds are only influenced by the dynamical variables, like position 
and velocity, of other birds in their immediate neighbourhood. Heterogeneity 
is introduced in the flock in terms of a depletion time after which a bird feels 
an increasing bias to move towards the landing surface. The stochastic nature 
of the bias ensures that the external and internal effects influencing the flock 
are also included. Through our model we have demonstrated a mechanism by 
which animals in large group can arrive at an egalitarian decision about the 
time of switching from one activity to another in the absence of a leader. Our 
results suggest that the coherence of the collective action of landing can be 
enhanced by the random perturbations. 
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